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CHAPTER I 


Introduction 

This work is focused on the development of a reliable 
fabrication technology for producing high quality 
Superconductor-Insulator-Superconductor (SIS) junctions for 
mixer applications in radio astronomy. Schottky barrier 
diode devices have generally been the standard in radio 
astronomy for low noise highly sensitive mixers in 
the millimeter and submillimeter wavelength regions. In the 
past couple of years, however, the lowest mixer noise 
temperatures, for frequencies below 115 GHz, have been 
obtained with SIS devices (see Fig. 1.1) [1], SIS mixers also 

require a local oscillator power level several orders of 
magnitude lower than conventional devices. The potential for 
conversion gain, predicted by John Tucker in his theory of 
quantum mixing [2], has been physically realized by several 
experimentalists [3,4] . Tucker's theory also predicts a 
minimum mixer noise temperature limited only by quantum noise 
(fiw/k). These advantages of SIS devices have prompted our 
research in this area. 

Chapter II presents a brief introduction to the 
phenomena of superconductivity, the concepts of "Cooper 
pairs", and quasiparticle excitations, and the development of 
the superconducting density-of-states. in Chapter III, basic 
SIS tunneling theory is developed and the typical SIS I-V 
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characteristics are examined. A theoretical modeling of the 
physical device parameters for use as optimal mixers is 
developed in Chapter IV. In Chapter V, initial fabrication 
results are discussed and a subsequent deterioration of 
device quality examined. The resolution of this degradation 
problem is presented in Chapter VI and potential avenues for 
future research are outlined in Chapter VII. 



CHAPTER II 


Superconductivity - The Phenomena 

Superconductivity is a term used to describe the quality 
of materials which exhibit a state change to zero resistivity 
when cooled below a critical temperature (T c ). The loss of 
resistivity is typically abrupt for pure samples, with T c 
dependent on the material. Investigation of the properties 
of metals at very low temperatures followed the discovery of 
helium liquif ication by Kammerlingh Onnes in 1908. 

The sharply contrasting zero resistivity of 
superconductors can be better understood by examining the 
specific heat C s of the transition. Figure 2.1 reveals a 
discontinuity in the specific heat at the critical 
temperature. The expression for the specific heat is given 
by 


C s = vT ds/dT, (2-1) 

where v is the volume per unit mass and s is the entropy. An 
entropy diagram Fig. 2.2 can therefore be obtained from Fig. 
2 . 1 . 

The specific heat can be separated into two parts - the 
crystal lattice component C x and the conduction electrons 
contribution C e ^. X-ray crystallographic investigations have 
shown that lattice properties are unchanged when 
superconductors are cooled below T c . The change in specific 
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S 



Figure 2.2 


Entropy for normal and superconducting 
states 



Specific heat for normal and superconducting 
states 
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heat between superconductors and normal states is therefore 
due to the conduction electrons. Since the entropy of the 
superconducting metal decreases with respect to the normal 
state, the electrons in the superconducting state must 
undergo a transition to a higher form of order. The 
electronic specific heat of metals in the superconducting 
state, far below T c , has been experimentally shown to have an 
exponential dependency on T given by 

C el = a exp (-b/kT) , (2-2) 

where a and b are constant. Such an exponential dependency 
is expected for materials with a gap in the energy spectrum of 
the electrons [5]. The existence of an energy gap is 
crucial to the theory of SIS devices and will be investigated 
further in this chapter. 

The phenomena of zero resistivity for constant direct 
currents indicates the absence of a voltage drop across the 
metal. For alternating currents an electric field exists 
such that resistivity is present. This a.c. resistivity, the 
general entropy diagram and many other characteristics of 
superconductors can be explained with a simple, physically 
appealing, "two-fluid" model. In this model the conduction 
electrons are separated into two species - "superelectrons", 
which induce the unusual superconducting properties, and 
"normal electrons", which behave like electrons in a normal 
metal. As the temperature is lowered below T_ the fraction 
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of superelectrons increases from zero until at absolute zero 
all are superelectrons. In the two fluid model, as the 
temperature is lowered below T c , some normal electrons 
undergo a transition to highly ordered superelectrons, 
explaining the lower entropy states for the superconducting 
material. 

Constant direct current must be carried by the 
superelectrons; however, for alternating current, the 
inertial mass of the superelectrons presents an inductive 
impedance, causing the superelectrons to lag behind 
the changing electric field. In an electric field the normal 
electrons are accelerated, presenting a resistive impedance. 
Bulk superconductors can therefore be represented by an 
inductance in parallel with a resistance. 

Further evidence of a gap in the energy spectrum of 
electrons is obtained from studies of reflection coefficients 
of the superconducting and normal states. A typical 
absorption curve (Fig. 2.3) [6] shows the existence of an 

absorption edge, for a superconductor, similar to the 
infrared radiation edge found with semiconductors. While the 
semiconductor infrared edge indicates an energy gap on the 
order of electron volts, the superconductor absorption curves 
imply a much smaller gap, on the order of milli electron 
volts. The two fluid model provides no insights for 
understanding this unusual phenomenon. 
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Another interesting phenomena that provides insight into 
the source of the energy gap is the characteristic of a 
"long-range order" by the superelectrons. It has been 
experimentally shown that the physical definition between 
superconducting and normal material cannot be sharp. The 
properties of the superelectrons can change significantly 
only over a distance greater than the "coherence length" - 
typically of the order of 10 -4 cm. This long range order 
implies that an unusual interaction between the electrons 
occurs. The fairly strong repulsive coulomb interaction can 
not create an energy gap. The mechanism must be attractive 
and much weaker to create the small energy gap measured for 
superconductors. In 1950, Frohlich proved theoretically that 
two electrons can be coupled together through an electron- 
phonon interaction. The electrons exchange momentum through 
the lattice in a virtual process, while the lattice is 
effectively undisturbed after the exchange. If the energy 
difference of the electrons, after the virtual phonon transfer, 
is less than the energy of the phonon, then the interaction 
will be attractive. This result, though proving the 
existence of a bound state between two electrons in a 
lattice, does not address a more realistic case of possible 
interaction between all of the conduction electrons. 

Before examining the BCS theory of superconductivity 
which followed Frohlich's work, a number of associated terms 
and notation must be developed. The simple case of the 
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electron distribution in a normal metal will first be 
examined, establishing the notation for the more complex 
superconductive case. 

The states that electrons can occupy ( (r)) are given by 
the solution to the Schrodinger equation, 

h 2 9 

- ^ V l(r) + V (r) tp(r) = e^(r), (2-3) 

where V(r) is the electric potential presented by the 
lattice. The solution to (2-3) can be given in the form of 
Bloch functions. With the well known f ree-electron model of 
metals, electrons can occupy Bloch states according to 

N n (e) = k/e (2-4) 

where N n (£) is the density-of-Bloch-states for a normal 
metal, k is a constant dependent on the electron mass and e is 
the energy of the Bloch state. N n (e) is shown in Fig. 2.4. 

The distribution of electrons with energy is described by the 
Fermi-Dirac statistics such that the probability of occupancy 
at energy e is given by 


' ' 1 + exp[ (e-e p )/k B T] 

where £p is the Fermi level. This function is shown in Fig. 
(2.5) for both the T = 0 ground state and the T > 0 cases. For 
T > 0, it is seen that there exists electron excitations for 
the case of a normal metal. 
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Density of Bloch status versus energy for 
a normal metal 



Fermi-Dirac distribution function 



Bloch state energies versus crystal 
momentum 
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The relationship between Bloch states energies and 
crystal momentum p, obtained with the f ree-electron model 
and the Schrodinger equation, is given by 



e 


and shown in Fig. (2.6). The e versus p diagram can be used 
to illustrate the phenomena of electron excitation from the 
ground state. Since the density of Bloch states presents a 
continuum of states in energy, any increase in temperature 
from T = 0 of the system, will result in electron 
excitations. It will be shown later in this chapter that 
this effect does not occur with superconductors, due to 
the energy gap that exists between the ground state and 
superconducting excitation states. Bloch states above and 
below e p will be designated by e e and respectively. With 
the system at T = 0, if conceptually one electron is excited 
from a Bloch state energy to e e ^, then the energy of the 

subsequent state above the ground state is called the 
excitation energy E. The excitation energy can be thought of 
as being composed of an electron excitation energy E e ^ and a 
hole excitation energy E^j. This definition of excitation 
energies is illustrated graphically in Fig. (2.7) with 
respect to Bloch state energies. For this conceptual case, 
the excitation energy is given by 


E " E ei + E hj = < e ei + e F> + < e F " e hj> 


(2-7) 
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For T > 0, the excitations are distributed continuously with 
energy, as is illustrated in Fig. (2.8). It is seen again that 
a continuum of excitation energies exists above the ground 
state energy E = 0. 

Returning now to the theory of electron-phonon 
interaction. Cooper, in 1956, built upon Frohlich's work 
by investigating the effects of adding two electrons to a 
metal at absolute zero. It is assumed artificially that the 
existing electrons in the metal obey 'normal' Fermi-Dirac 
statistics, which presents a step-function at absolute zero (Fig. 
2.5). If there is an electron phonon interaction, the two 
electrons can be thought of as scattering each other such 
that their individual momenta are continually changing with 
time. At any instant of time, the two electrons can be 
viewed as having some momentum p i and pj. The continual 
scattering of the two electrons can be thought of as being 
composed of a collection of such scattering cases (or 
"events") under the constraint that the total momentum is 
constant. Cooper showed that potential energy is minimized 
for the system with the greatest number of permitted 
scattering combinations, i.e. the total number of possible 
(p if pj) combinations. If the electrons with equal but 
opposite momentum are paired, then the maximum number of 
allowed scattering pair states are available, such that the 
system energy is minimized. The pairing of two electrons, p.l 
and -Pil (where "1" denotes spin up) is now known as a Cooper 



( 



Excitation energies defined in terms of 
Bloch state energies 
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pair. Cooper showed that this type of pairing yields the 
largest decrease of potential energy. If two electrons, which 
previously were in normal states made a transition to such a 
scattering state, then the decreased potential energy would 
offset the increased kinetic energy such that the total 
energy of the system would decrease. 

A single electron can be represented by a quantum 
mechanical waveform r ,p) , where ip describes the 

probability that the electron has momentum p and is located 
at position r. The probability function for two uncorrelated 
electrons with momenta p^ and pj and position r is described 
by the wavefunction 

<l> ' (t/P if Pj) = ^ ( r ,Pi)^ (r,Pj) (2-8) 

The Cooper pair state, which consists of two paired electrons 
in a mixture of wavef unctions with equal and opposite 
momenta, is described by the wavefunction where 

-Mr) = l a i j4>(r,p i l,-p i l) , (2-9) 

ij 

If the system energy is lowered by the formation of one 
pair, then it would seem logical for many electrons to make 
the transition. The probability of finding a given pair 
state (Pj +,-pj + ) empty at any given time decreases as more 
Cooper pairs are formed with p > pf Open momentum states 
become increasingly unavailable to the Cooper pairs for 
occupation. Since reduction of system energy is a function 
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of the number of scattering processes available (i.e. the 
number of possible (Piz-p^) states available for the next 
scattering event per electron pair), there exists a tradeoff 
between increased kinetic energy and reduced potential energy. 
There is therefore a limit to system energy reduction by 
forming Cooper pairs with p > p f . 

In 1957 Bardeen, Cooper and Schrieffer developed a 
theory of superconductivity which addressed this question of 
many pair interaction [7]. The crucial assumption of 
the"BCS" theory is that pairs influence each other only by 
making unavailable, for the next scattering event of the other 
pairs, the momentum states which they are occupying. In the 
BCS theory, at absolute zero all the conduction electrons 
form Cooper pairs with the wavefunction $ given by 
eqn. (2-9). These pairs are described by Bose-Einstein 
statistics and exist at one energy value in the same quantum 
state. Some of these pairs form with p > p f as described 
earlier such that the system energy is lowered. The 
probability that a two electron state (p^l, -p^l 

is occupied, is given by 



£ i " £ F 


((e.-e F ) 2 + A 2 ) 172 


( 2 . 10 ) 


where 


(Pi 2 /2m), and A is given quantitatively by 


A (T) ^ 

A (0) T 


A (T) 


(0) tanh 


( 2 - 11 ) 
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Figure 2.9 shows the probability h^ that the double-electron 
state is occupied as a function of momentum at T = 0. The 
superconducting state has the unusual property of possessing 
vacancies below p^ and occupied states above Pf for the T = 0 
case, unlike a normal metal. 

From BCS theory, the energy necessary to break apart a 
Cooper pair is given by 

2 2 2 2 

E b = ((e i -e F ) 2 + ) + ( (e h -e p ) + /T) (2-12) 

where e e and are Bloch state energies above and below £p 
respectively, as defined earlier. The theory therefore 
predicts a gap of 2A in the excitation energy of the 
superconductor, centered about the fermi energy. This energy 
gap, as noted above, was seen experimentally in specific heat 
and absorption studies but could not be explained with the 
simple two-fluid model. The separation of a Cooper pair 

+ implies that the electrons no longer form that 

pair with equal and opposite momentum. If an excitation to 
momentum state p^ + occurs, then its Cooper pair partner is 
absent from state -p^ . These excitations (called 
"quasiparticles") no longer participate in the scattering 
interactions - they act almost like free electrons. The 
splitting of the Cooper pair increases the potential energy 
of the system by its binding energy. With an excitation in 
momentum state p^i but the state -p^l empty, the pair state 
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Figure 2.9 

Probability that the double-electron state 
is occupied 



Superconductor energy gap dependence on 
temperature 
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(Pi+,-Pi+) is unavailable to the other pairs for occupation 
in the next scattering event. The energy of each electron 
pair in the total system is therefore raised. These two 
energetic effects cause the noted gap in the excitation 
spectrum. As the temperature is raised and quasiparticles 
are created, the decrease in scattering interactions yields a 
decrease in the energy gap as shown in Fig. 2-10. The gap 
energy continues to decrease as the temperature is raised 
until the critical temperature is reached and A = 0. It is 
important to note that A does not substantially decrease 
until the temperature is raised above T c /2. 

The excited electrons are named quasiparticles because 
of the unusual properties of these states. Consider again 
Fig. 2.11 which shows h^ 2 , the probability that the double- 
electron state is occupied. If a quasiparticle is excited 
with momentum p^, then the state -p^ is empty. If p^ is much 
greater than pf, then the excitation exists in the state p^, 
which previously had a very small probability for occupancy 
(hi 2 << 1). The electron therefore exists in a state which 
was previously unoccupied, such that the quasiparticle acts 
like an electron (and is called a quasielectron). Similarly, 
if Pj is much less than p f , then the excitation exists in the 
state pjf which previously had a very large probability for 
occupancy (h^ >> i). The electron therefore exists in a 
state which was previously occupied. The corresponding state 
-p j , which previously was occupied is now unoccupied, such 
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Probability of pair occupancy in the 
superconducting state 



Figure 2.12 

Permitted excitation energies for a 
superconductor with energy gap A 
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that the excitation acts like a hole (and is called a 
quasihole). If is near pj, such that h^ is a fraction of 
unity, then the quasiparticle cannot be regarded totally like 
an electron or a hole. Instead the quasiparticle acts as a 
mixture of both. 

The permitted excitation energies (E) versus crystal 
momentum, as given by BCS theory, are shown in Fig. 2.12 for 
T = 0 and T > T c . Away from p f , the E versus P diagrams for 
the superconductor and the normal case are equivalent. Note 
the gap in the excitation energy of energy A for the T = 0 
case. In contrast to the normal metal case, excitations can 
only occur in the superconductor case for energies larger than 
A. 

Earlier in this chapter, the density-of-Bloch-states for 
a normal metal was given. In order to analyze the 
tunneling phenomena of quasiparticles in a SIS junction, the 
density-of-excitation-states must also be known. For the 
case of a normal metal, it is evident from Figs. 2.7 and 2.8 
that, for every value of p, there exists a corresponding 
Bloch energy and a corresponding excitation energy. For the 
Bloch states, in an incremental interval of momentum dp, 
there is a corresponding range of Bloch state energies de. 

For the excitation states, in this same interval of momentum 
dp, there is a corresponding range of excitation energies dE. 
These two cases are shown in Fig. 2.13. Since the number of 
Bloch states and excitation states in this interval dp, are 
equal, the number of states in de equals the number of states 
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in dE. The density-of-normal-excitation-states, N n (E) , is 
therefore given by 

N n (E)dE = N n (e)de. (2-13) 

For the normal metal, the excitation energy is just 
E = |e - e p | , as discussed above, such that dE = de and the 
density of excitations for a normal metal, N n (E), is simply 
equal to N n (e), the density-of-Bloch-states. Though the 
approach for deriving the density-of-excitation-states- 
function of the superconductor is simlar, the resulting 
function will be shown to be quite different. 

The derivation of the density-of-excitation-states for a 
superconductor will now be examined. From Fig. 2.8, it is 
evident for the superconductor that for every value of p 
there exists an excitation state such that (2-13) is modified 
to 

N s (E)dE = N n (e)d . (2-14) 

Because there exists a gap about the Fermi level in the 
excitation energy, the simple relationship de/dE = 1 does not 
hold for the superconductor. In fact, the value of the 
density-of-superconducting-excitation-states is expected to 
be equal to zero for E < A. However, far above or below 
p f , the E versus p diagrams for the superconductor and the 
normal case are equivalent such that the superconducting 
density-of-states N s (e) is identical to that of a normal 
metal N n (E). Near the energy gap of the superconductor, 
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Graphical analysis of N S (E) 


Figure 2.15 
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however, N S (E) will be quite different from N n (E). 

The form of N S (E) will now be illustrated graphically. 

In Fig. 2.14, an expanded view of the superconductor 
excitation energy near Pf is shown. Since the energy gap has 
such a relatively small value, when we expand the energy axis 
to show A graphically, the E versus p curve away from pj 
becomes nearly vertical. In order to illustrate the origin 
of the unusual superconducting density-of-states, the scale 
of the E versus p diagram for Fig. 2.15a is distorted. 

Figures 2.15b - 2.15d show e versus p (normal case), N g (E) 
and N n (e) respectively. For momentum values far from 
Pf, an incremental interval dp^ is shown to map into a range 
of normal Bloch state energies (de n ^) and superconductive 
excitation energies (dE s ^). Since N n (e)de nl = N s (E)dE sl , the 
density-of-excitation-states for the superconductor can be 
constructed graphically, as N n (E) is shown. From Fig. 2.15d, 
the shaded area, under the known density-of-states curve of 
the normal metal for the de nl interval, represents the number 
of normal states for this interval and is denoted N n ( £ )de nl . 
In Fig. 2.15c, N S (E) (corresponding to the dE gl mapping) is 
constructed such that the shaded area N g (E)dE s ^ is equal to 
N n (e)de nl . Since the chosen momentum interval dp-^ is not 
close to pj, the density-of-excitation-states construction 
for the superconductor is similar to the normal case as seen 
in Figs. 2.15c and 2.15d. If another momentum interval (dp 2 ) 
closer to p f is chosen, the energy gap will strongly alter 
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the density-of-states function for the superconductor. Due 
to the decreasing slope of the E versus p superconductor 
curve, dp 2 intersects a superconductor excitation interval 
dE S 2 which is much smaller than de^. In order to construct 
a N g (E)dE S 2 e< 3 ua l to N n (£)de n 2 , the superconductor density- 
of-states must increase greatly with respect to the normal 
function. As p approaches p^, the interval of excitation 
energies decreases to zero such that the density-of-states 
function N g (E) approaches infinity. Note that no excitation 
states exist in the energy gap; they can be thought of as 
having been pushed above E = A. N g (E) , for E < E p , can be 
graphically derived in a similar manner. From BSC theory, 
the excitation density-of-states for a superconductor is 
given by 


v e > 


N (E) E 


(.2-15) 


If the undistorted energy scale of Fig. 2.14 is used, 
then the normal density-of-states curve away from p^ is seen 
to be nearly constant for tens of energy gaps. For this work 
the energy range of interest is less than + 10 A. N n (E) can 
therefore be approximated by N n (0), such that (2-15) 
becomes 


V E) 


N (0) E 
<E 2 -A 2 > 1/2 


(2-16) 
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N S (E) is shown to scale in Fig. 2.16. In order to show N S (E) 
graphically for this work's purposes, the scale will again be 
distorted as shown in Fig. 2.17. Here the density-of-states 
is shown in what is known as the "semiconductor" 
representation, where the excitation energy is positive for 
both directions from the origin. The top half denotes 
quasielectron excitation energy and the lower half, quasihole 
excitation energy. 




CHAPTER III 


SIS Tunneling Theory 

As an introduction to SIS tunneling theory, the case of 
normal metals (NIN) will first be examined. For an insulator 
which presents a relatively larger barrier, the thermal 
component of current over the potential barrier is negligible. 
Current transport is due to tunneling through the oxide with a 
tunneling probability P which can be given by a quantum 
mechanical analysis of the particular barrier shape. In 
general, however, P decreases exponentially with barrier 
thickness and increases linearly with area. The actual 
characterization of a barrier is often limited to a tunneling 
resistance which contains P implicitly. 

Figure 3.1a is an electron population diagram for a NIN 
junction at T = 0. As discussed above, since the energy range 
of interest is on the order of millielectron volts, N^(£) has 
been taken equal to N n (0). For T > 0 the electron occupation 
is therefore described by the Fermi-Dirac distribution f(E), 
such that the resulting population is given by Fig. 3.1b. 
Tunneling transitions must of course be made horizontally on 
these diagrams to conserve energy. The current due to 
tunneling from left to right in an energy interval dE, is 
found by simple analysis to be proportional to 


Il+2 aP[N 1 (E-eV)f (E-eV)] [N 2 (E) (1-f (E))]dE 


(3-1) 
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while the current from right to left is proportional to 


1 2 + 3.01 P [N x (E-eV) (1-f (E-eV) ] [N 2 (E) f (E) ] dE. (3-2) 

The total current is found by integrating over all energies. 
For the NIN case, the current expression is 


or 


I = BP N 1 (0)N 2 (0) / (f (E-eV)-f (E) )dE (3-3) 

all E 


I = BP N 1 (0)N 2 (0)eV (3-4) 

where B is a constant which is related to the geometry and the 
integral has been approximated by eV. Equation (3-4) can be 
written as 


I = Gnn V, (3-5) 

where Gnn is the conductance of the ohmic device in the normal 
state. 

Figure 3.2 shows an excitation tunneling diagram for a 
NIN junction. This form of tunneling diagram will be very 
useful in understanding the physics of superconductive devices 
later in this chapter. The NIN case will serve as a simple 
introduction to the concepts involved. In Fig. 3.2, a 
negative bias has been applied to the left side (metal one) 
such that the zero excitation level of metal one is raised by 
an energy eV relative to metal two. The tunneling of an 
electron from metal one to metal two is represented by a hole 


i eV-E-j+E2 
i r 


Y. 

P« 


P 


Figure 3.2 Excitation tunneling diagram for a NIN 

junction 



Figure 3.3a Excitation tunneling diagram for a NIS 

junction 



Graphical demonstration that U- 2 +Uu 2 = 


Figure 3.3b 
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excitation (p < p f ) of energy E^ in metal one and an electron 
excitation (p > p f ) of energy E 2 in metal two. The applied 
bias results in excitation energies such that E i + E 2 = eV. 

The case of a SIN excitation tunneling diagram is shown 
in Fig. 3.3a, where an electron of energy Ej can tunnel into 
state Pj or p k , both states at energy E 2 . For the T = 0 case, 
the probability that state p k is empty is given by 
u k = 1 -h k . Since h z ^ has symmetry about p f it is also 
evident that h 2 j = 1 - h 2 k or u 2 k = h 2 j (Fig. 3. 3b). Since the 
sum of the probabilities for occupancy and non occupancy of a 
given state is unity (h 2 j + u 2 j = 1) it follows that 

Uj 2 + u k 2 = 1 (3-6) 

Equation (3-6) reveals that the total probability of a 
vacant state is unity. Though the electron may actually 
tunnel into either state Pj or p k , the current can be 
conveniently modeled by describing all states below p p as full 
and those above p F as vacant. In this scheme, the electron 
tunnels into the vacant state p k . The semiconductor 
representation for the electron population diagram of a SIN 
junction at T = 0 is therefore shown in Fig. 3.4. In this 
representation, electrons tunnel from the normal metal only 
into states above the energy gap, though physically states 
both above and below Pf are involved. 

Figure 3.5 shows an excitation tunneling diagram for a SIS 
junction. The representation of states in the superconductor 
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as empty above p F is used for both superconductors. Ideally, 
there are two tunneling processes. 

In the first case (for T = 0) a Cooper pair may separate 
into a hole excitation and an electron which subsequently 
tunnels through the barrier into the other superconductor as 
an electron excitation. The minimum applied bias for 
tunneling occurs when the bias is equal to the sum of the two 
gap energies. This process is shown in Fig. 3.6 (for T = 0) 
using the semiconductor representation. The second process is 
the tunneling of an existing excited electron across the 
junction, occurring for T > 0. The general semiconductor 
representation of a SIS junction for T > 0 is shown in Fig. 

3.7. BCS theory shows that the distribution of single 
electron states in a superconductor is of the Fermi-Dirac form 



(3-7) 


where E D is the excitation energy corresponding to momentum 

r 

state p from Fig. 2.12. The excitations above and below the 
gap energies in Fig. 3.7 (corresponding to quasiparticle 
occupation of states about p F ) are therefore distributed 
according to the Fermi function. The shaded areas (for E e ) can be 
thought of as potential tunneling electrons and the unshaded 
areas as open states for tunneling electrons. Alternatively, 
the shaded areas (for E^) are empty hole states and the unshaded are 
are empty electron states. 
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Consider first the T = 0 case for two superconductors with 
gap energies A1 and A2. When a positive bias is applied to 
metal two, electron tunneling is not permitted until 
eV = A1 + A2, as shown graphically in Fig. 3.8. At this 
potential, the singularities in both density of state 
functions are aligned such that (with our semiconductor 
model) the tunneling of electrons from metal one to metal two 
increases substantially. Since the superconductor density- 
of-states curve approaches the normal case for energies 
larger than A, the current approaches the normal curve for 
large biases. The T = 0 I-V curve, for this SIS case, is 
shown in Fig. 3.9. 

The general current expression, analagous to (3-3) is 

I = 8P /N sl (E-eV)N s2 (E) {f (E-eV)-f (E) }dE. (3-8) 

If eqn. (2-16) is substituted, the expression simplifies to 

I = N sl (E-eV)N s 2 (E-eV) {f (E-eV)-f (E) }dE (3-9) 

where N' SX (E) is defined as 

N sx (E) = 2 E T 1/2 (3_10) 

SX (E 2 - A t) ' 

Equation (3-9) can be solved directly when T = 0 to give 

Fig. 3.9; however, the T > 0 case can only be solved numerically. 

For the T > 0 case, excited quasiparticle tunneling is 
permitted through the barrier for biases less than 
Ai + a 2 * When a negative potential is applied to metal one, 
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the subgap current will increase until eV = A 2 - A-^. At this 
bias, the upper edges of both energy gaps are aligned as shown 
in Fig. 3.10. With the semiconductor model, for larger 
biases, the thermally excited electrons above the gap and 
holes below the gap tunnel into a decreasingly available 
density-of-states. The current, therefore, decreases until 
eV = a ^ + A 2 , yiel<3i n 9 a negative resistance region. With a 
SIS homojunction, where A x = A 2 , this phenomena does not 
occur, such that this negative resistance region is peculiar 
to heterojunctions. For eV = A 1 + A 2 , the tunneling current, 
as in the T = 0 case, increases substantially and approaches 
the normal curve for larger biases. A typical T > 0 ideal I-V 
curve is shown in Fig. 3.11. 

Though the exact solution of the current equation (3-9) 
requires a numerical integration, a number of theoretical 
properties of the I-V structure have been derived from BCS 
theory. The quasiparticle current discontinuity at 
eV = A1 + A 2 is given by 

Al ss = [A 1 A 2 l 1/2 [l-exp[-(A 1 + A 2 )/kT]]/ 

(3-11) 

[1 + exp(-A 1 /kT) ] [1 + exp(-A 2 /kT)] 

For temperatures below T c /2, such that the energy gap has 
nearly obtained a A(0) value, AI gs can be approximated by 

AI SS - SjHF [Ai 4 2 ] 1/2 (3-12) 
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In addition to quasiparticle tunneling, pair tunneling can 
occur in a SIS junction. Pair tunneling can be described by 
the Josephson equations 


and 


I = I c sin <P 


d<f> _ ? eV 
dt fi 


(3-13) 


(3-14) 


where <f> is the phase difference across the junction of the two 
pair wavefunctions. Two separate phenomena, the dc and ac 
Josephson effect, are described by these equations. The ac 
effect occurs when <J> changes at a constant rate in proportion 
to the applied voltage. This results in an alternating 
current at the frequency 2 ^ V^ c . The dc effect occurs when 
$ is constant such that pairs can tunnel through the barrier 
even though the voltage across the junction is zero. When the 
critical current (I c ) is exceeded, a voltage develops across 
the junction and <p begins to change. 

The dc Josephson critical current is often suppressed by 
stray magnetic fields and unwanted noise. The low temperature 
form of I c , derived from BCS theory, is 


T = — - — — — — 
c eR A , +A _ 
n 1 z 

With eqn. (3-12), the ratio of the critical current to 
quasiparticle current jump is found to be 


(3-15) 



2/a7a7 


(3-16) 
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such that the ratio is theoretically less than one, except for 
the case. If Ic is suppressed, the value can 

theoretically be obtained from Al ss and eqn. (3-16). 

Figure 3.12 shows a theoretical I-V curve for two 
electrodes (A(Nb) = 1.44 mV and A(PbBi) = 1.7 mV ) which are 
used in this work. With these materials, the quasiparticle 
current jump AI SS = 49.15 yA is theoretically only slightly 
larger than the critical current I c = 48.98 yA, for R n = 50 A . 
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Theoretical I-V curve for an SIS junction 
A-^= 1.44 meV, ^ = 1.7 meV 


Figure 3.12 


CHAPTER IV 


Theoretical Modeling 

This work is focused on developing an optimal fabrication 
technology to produce high quality SIS junctions for mixer 
applications. In this chapter, the mixing requirements on 
junction parameters are stated and the theoretical modeling of 
these parameters is subsequently developed. The pair of 
electrodes utilized in this work (Nb/PbBi) are examined. 

Ideally, the I-V curve of an SIS mixer should look like 
that of a switch - the subgap current below V = + A2 should 

be small and the quasiparticle current jump should be abrupt. 
Unfortunately, the ideal SIS I-V characteristics (Fig. 3.11) 
are usually altered by numerous nonideal effects. For real 
devices, the subgap current is found to be higher and the 
current jump less abrupt than theoretically predicted. These 
effects will be discussed later in this chapter. 

The abrupt nonlinearity of the SIS quasiparticle current 
causes an unusual effect in the mixing process. While the 
D.C. nonlinearity ideally occurs over a small percentage of a 
millivolt, the voltage equivalent of the photon energy at 115 
GHz is 0.475 mV. Unlike the case of Schottky mixers at this 
frequency, an application of L.O. power will "significantly" 
add to the D.C. bias due to the phenomena of photon assisted 
tunnelling. The effective excitation energies are modulated 
by the applied L.O. field. This phenomena is illustrated 
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(for the T=0 case) in Fig. 4.1b with a semiconductor diagram. 

In this diagram three cases are illustrated: no L.O. field, 

photon absorption and photon emission. With no application 
of a L.O. field, it is seen that an applied bias of 
(A 2 + A 2 )/e is necessary for current flow. For the second 
case, however, the excitation energy of electrons, which have 

t 

absorbed photons of frequency V L0' has increased by hu> L Q. The 
associated energy levels in the semiconductor diagram have 
been shifted by ho). It is seen that a smaller applied bias of 
(A i + A 2 - ■fito)/e is therefore required for current flow. For 
the third case, the excitation energies of electrons, which 
have emitted photon of frequency V LQ , have decreased by 
'fiaiLo- It is seen that a larger bias of (A-^ +A 2 + ‘'Kw)/e is 
therefore required for current flow. 

The pumped I-V curve can therefore be thought of as a 
superposition of unpumped I-V curves, each one shifted on 
the voltage axis by an integer number of photon step 
increments (of voltage liw/e) . A typical SIS I-V curve with 
and without L.O. bias is shown in Fig. 4.1b. 

If the instantaneous voltage is driven below a certain 
voltage (called the dropback voltage, V^) the voltage across 
the junction drops to zero. The critical current 
of the junction must then be surpassed before the device can 
return to its quasiparticle branch. 

This Josephson effect noise must be avoided in SIS 
mixers; i.e., the D.C. bias should exceed V^ + V L0 . However, 




representation of photon emission 
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the optimum D.C. bias is approximately V^ c ~ Vg -^fiu)/2e while 
the optimum L.O. voltage is V LQ = -fiuj a /e (where a should be 
relatively constant) [9]. The condition for avoiding this 
Josephson effect noise, therefore, becomes 


V 

gap 



(a + |) 


(4-1) 


With a Josephson effect analysis, the dropback voltage is 
found to be proportional to [10] 


v d * 


(wR n c ) 


1/2 


(4-2) 


In contrast to the case of Schottky diodes, a RC cutoff 
frequency does not exist for the superconductive device. 
Capacitance is in fact desirable for SIS mixer operation in 
order to provide a lossless termination for the undesirable 
upper sidebands. If the device capacitance is too large, 
however, the signal bandwidth is reduced and optimal R.F. 
tuning is complicated. A value of Q=4 (Q ~ wR n C) has been 
demonstrated by Feldman [11] and Kerr [12] to be an effective 
compromise where a normal resistance of R n * 50 ft is typically 
chosen as a compromise between RF and IF matching 
requirements . 

Since Q ~ 4, from eqn. (4-2) it is seen that the dropback 
voltage is proportional to ojl/2. There is, therefore, a 
maximum frequency of operation for which the inequality of 
eqn. (4-1) is obeyed. The upper frequency limit of an SIS 
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mixer with Vg ~ 3 mV has been numerically estimated to be 
[13] 300 GHz. 

This upper frequency estimation is not/ however/ an 
absolute limit. A number of different techniques exist for 
supressing the Josephson critical current and reducing V^. 
[14/15] The potential for the use of Josephson effect mixers 
[16] at higher frequencies has also been largely unexplored. 

For a given set of electrode and barrier materials/ there 
are two physical variables to be chosen - barrier width and 
area. The junction capacitance and normal state resistance R n 
are set by a choice of these variables. A value of Q = 4 
results in a capacitance 


c 


6.366 x 10 5 
f R 

n 


fF/ (GHz - ft) . 


(4-3) 


For R n = 50 ft and f = 115 GHz f the desired junction 
capacitance is 111 ff. 

Experimental optimization of a given SIS device design 
requires a method for obtaining a value for the R n C product. 
Experimentally, R n is relatively easy to determine, however, 
the measurement of capacitance can be quite involved and 
difficult. If the tunneling barrier thickness, area and 
dielectric constant are known, then the desired capacitance 
can be calculated from 


C 


d 


(4-4) 
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where e r is the barrier dielectric constant. However, the 
barrier thickness is very difficult to measure. A value for 
C, however, can be obtained indirectly from the critical 
current. Since R n increases exponentially with d, and with 
eqn . (3-15), J c can be written as 

J c = J Q e" Kd (4-5) 


where K is a material tunneling constant. With eqn. (4-4), 
eqn. (4-5) can be rewritten as 



InJ 


'Ke e ' 
r r 




lnJ c 


(4-6) 


where C A is the capacitance per unit area of the junction. 
Experimental results for Nb/PbBi junctions [17] have been used 
to fit C A to within approximately 10% for 100 < J c < 

5,000 A/cm^) yielding 



10.5 ^ - 0.45 

pf pf 


In (- 


A/cm^ 


( 4 - 7 ) 


where J c is the experimentally observed value of the critical 
current density. With eqn. (4-7), a value for C can be found 
with a measurement of I c and A, such that a ujR n C product is, 
therefore, obtained. 

The interrelationship between Q,R n ,J c ,C,A and d will now be 
examined in detail. With eqn. (4-7), K and J 0 can be 
identified such that solving eqn. (4-5) for d yields 
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d 


J 

23.3 - In ( — 2__) 
A-crtt 

8.43 x 10 3 


ym 


(4-8) 


J c as a function of d, for Nb/PbBi junctions is shown in Fig. 
4.2 where J c is seen to depend strongly on d. In order to 
fabricate devices with a desired J c value, a tight control on 
the barrier thickness is required. 

It has been shown that the expected BCS value of the 
critical current must be multiplied by an empirical 
correction factor (0.64 for Nb/PbBi [18]) to correspond with 
experimental results. For the experimental measurement of J_, 
(3-15) is, therefore, modified to 


J 


c 


0.64 tt A 1 A 2 

AeR A„ 

n 12 


(4-9) 


If eqn. (4-9) is multiplied with eqn. (4-4), then Q can be 
written as 


Q 


(jOtt (0.64) 
J 

c 


A 1 A 2 
a i+a 2 d 


£ £ 

r o 


(4-10) 


With eqn. (4-8), it follows that Q and J c are inversely 
related. (J c versus Q is shown in Fig. 4.3). For a choice of 
j c (or Q) C can be varied such that from eqn. (4-7) values of 
A can be calculated and values of R n finally obtained from 
eqn. (4-9). With this method, the range of values of the 
parameters R n ,C,A and J c (or Q) , yielding desirable RF 
characteristics, can be found. In Fig. 4.4, lines of constant 
Q and A are shown on a R n versus C plot for a Nb/PbBi junction. 
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Figure 4.3 

J c versus Q, 

(Q = w R n C) 




57 


Acceptable junction characteristics fall between the Q = 3 and 
Q = 5 curves for R r between 30 and 70 ohms. Tabulated values 
of R nf C,A and J c have been listed in Table 4-1. It is seen 
for a given A, that C is almost constant, decreasing slightly 
for increasing R and Q. This is understood by noting that a 
very small increase in barrier thickness will result in a much 
larger increase of R n and decrease of J c values, yet the 
junction capacitance will decrease only slightly. Further 
insight can be obtained from Fig. 4.5, where constant R n 
curves of 30 and 70 ohms are plotted versus A and C. Again it 
is seen for a given area, that the capacitance is a weak 
function of normal resistance. 

Control of the fabrication process, as has been mentioned 
above, is of critical importance if the desired Q values are 
to be achieved. In Fig. 4.6, constant lines of R n are shown 
on a Q versus A plot where acceptable characteristics fall 
within the shaded area). The desired values of junction area, 
therefore, range from 0.4 to 1.6 square microns. Only a small 
range of areas, from 0.7 to 0.9 pm 2 , permit fabrication of 
devices with the full range of desired Q. 

Definition of areas less than 1 um 2 is difficult to obtain 
utilizing a planar geometry with standard photolithographic 
techniques. Special fabrication techniques such as shadow 
masking or an edge geometry are needed to achieve smaller 
geometries without resorting to electron beam lithography. 
Mixer use at higher frequencies necessitates even smaller 
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A(»m) 



Figure 4.5 


Lines of constant R n are shown on an A 
versus C plot 
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geometries. With Q = 5 and R n = 50 ft, a device designed for 
operation at 220 GHz and 350 GHz requires device areas of 0.48 
and 0.29 square microns respectively. 

An alternative to small area junction fabrication, is to 
utilize series arrays of junctions. If identical junctions 
are used, each junction with an area equal to N A, then the 
array presents a Q which is equivalent to that of a single 
junction with area A. Junctions with larger areas, which are 
easier to fabricate, can, therefore, be used to obtain the 
critical Q value. Since voltage and current scale with N, the 
use of arrays increases the saturation level of the device. 
Though arrays present potential advantages, for optimum use 
the array dimension should be small compared with A./4 and all 
junctions must have approximately the same R n [19] • These 
conditions, especially the uniformity of device 
characteristics, can prove difficult to achieve. 

In addition to the parameters discussed previously 
(Q,R n fA) the degree of deviation from the ideal I-V 
characteristics is very important for optimum mixer quality. 
The expected low subgap current and abrupt quasiparticle 
current jump are typically not fully obtained with actual 
devices. Often a ratio of the small signal resistance below 
the gap voltage to the small signal resistance above the gap 
voltage is given as a junction quality parameter. For 
Nb/Pb-alloy junctions, the chosen ratio is R(2mv)/R(4mv) or 
r2,4, such that the values collar the turn on voltage by 
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approximately lmv. For the Nb/PbBi case (Fig. 3.12) the ideal 
value of R2,4 (at 4.2K) is approximately 55. [20] 

There are three critical physical components in realizing 
a high quality SIS junction - the base electrode, insulating 
barrier and counterelectrode. The crucial step in device 

o 

fabrication is the formation of the regions between 100 A on 
either side of the barrier. A high quality base electrode 
material can be easily disturbed by the chosen cleaning 
process used prior to barrier formation or by formation of the 
barrier itself. The surface damage may reduce the T c and the 
gap energy or even change the material to a normal metal, 
while the effects can occur locally or over the whole junction 
area. The reduction of gap energy and formation of states 
within the gap will result in a distribution of effective turn- 
on voltages such that the theoretically abrupt quasiparticle 
current jump can become quite rounded and render the device 
useless for mixer purposes. 

Figures 4.7 a and b schematically show the effects of two 
cases of surface damage to the base electrode, assuming an 
ideal barrier and counterelectrode. In Fig. 4.7a, the damage 
has occurred in varying degrees across the surface, resulting 
in locally fluctuating values of gap energies. This device 
can be thought of as being composed of a number of junctions 
in parallel, each with a different Ai + Aj turn on voltage. 

The resulting I-V curve will have a low subgap current for 
voltages below the smallest effective gap energy. Above this 
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Figure 4.7a 


SIS I-V curve, for a case of a graded 
distribution of effective turn on voltage 



Figure 4.7b 


I-V curve for a NIS junction 
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voltage, however, the once abrupt turn on voltage will become 
rounded, due to the graded distribution of effective turn on 
voltages. In Fig. 4.7b, the damage is so severe that the 
material no longer has an associated gap energy near that of 
the barrier, such that the junction could be modeled as an NIS 
device. 

Unfortunately, damage to the base electrode often affects 
the subsequently deposited barrier material, yielding a leaky 
tunneling barrier. Numerous materials and methods exist for 
barrier fabrication - the most suitable selection depends on 
the choice of base and counterelectrode materials in addition 
to the device geometry and method of fabrication. The optimal 
barrier is of uniform thickness and absent of potential 
tunneling states. A leaky barrier typically results in an 
increased current (subgap current) for all voltages below the 
gap energy. 

Damage to the counterelectrode can result in detrimental 
effects mentioned previously for the base electrode. Damage 
to the barrier can also occur depending on the method of 
counterlectrode deposition and composition of the barrier. In 
the case of damaged superconducting materials, with a 
heterojunction the suspected material can be identified by 
examining the A1 - A 2 cusp and the A1 + A 2 turn on voltage. 

If the subgap current is fairly high, or the turn on voltage 
too rounded, accurate identification of the sum and difference 
of the gap voltages can be difficult. 
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The device fabrication techniques chosen to accomplish all 
these steps are of extreme importance; the complexity of this 
challenge will be considered in the following chapters for 
Nb/PbBi junctions. 



CHAPTER V 


Initial Fabrication Results 

The initial work on this project was undertaken by J. 
Upshur and Dr. R. Mattauch. An edge geometry was eventually 
chosen, utilizing the submicron thickness available with a 
thin Nb film and the less than 2 photolithography 
capability of the Semiconductor Device Lab to achieve a sub 
square micron junction area. A considerable amount of 
attention by other researchers, however, had focused on lead 
or lead alloy junctions [21,22,23]. One of the material 
advantages of Pb over Nb is that the dielectric constant for 
its native oxide is approximately 3-4 times smaller than that 
of The desired Q of 4 can, therefore, be obtained 

with a larger junction area for the Pb technology, obviating 
the need of an edge geometry. A planar geometry is typically 
chosen with a barrier formed by thermal oxidation in oxygen 
or even in the atmosphere. 

Though relatively promising results below 115 GHz had 
been obtained with Pb alloy mixers, [24,25] several 
disadvantages are associated with this technology. It is 
well known that pure lead films degrade mechanically with 
thermal cycling. The existing mechanical strain often 
results in damage to the oxide tunnel barrier, or even an 
electrical short between the two superconducting electrodes 
[26]. Though it has been shown that some Pb alloys (deposited 
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by observing precise schedules) have improved thermal cycling 
characteristics with respect to mechanical 
strain, Pb alloy devices have generally demonstrated a 
gradual drift in device characteristics with time [27]. Pb 
alloy films are difficult to pattern for small line widths. 

In addition, the films are typically not mechanically durable 
enough to sustain a simple photolithographic step. The use 
of the overlap technology associated with the Pb junctions, 
for frequencies greater than 115 GHz, would also prove very 
difficult because of this size limitation. 

The choice of Nb as a base electrode was made for 
several important reasons. Nb is a mechanically rugged 
refractory metal which is physically compatible with 
photolithography and numerous patterning steps. Nb devices 
have also demonstrated superior mechanical and electrical 
stability properties to lead devices. 

Though the edge geometry, as initially planned, utilizes 
both Nb and a Pb alloy, it was thought that this structure 
would be a promising step toward an all Nb junction. A 
schematic drawing of the Nb/Pb alloy device is shown in Fig. 
5.1. The technology utilized draws from earlier work at IBM 
and Cornell University in addition to considerable 
development at the University of Virginia. A general summary 
outline of the fabrication technology developed until January 
1982 is presented below. 


QUARTZ 


Figure 5.1 Schematic of a Nb/Pb alloy edge contact 

device 
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The substrates used are fused Suprasil II quartz 
(Amersil, Inc.), 19 mm x 19 mm x 7 mils thick, chosen for 
their superior R.F. properties. The Nb base electrode film 
was deposited by RF sputtering in an argon plasma at the 
University of Virginia Physics Department. The Nb target and 
substrate were sputter etched clean prior to deposition of 
the approximately 2,000 A thick film. The system base 
pressure was typically 2 x 10~^ torr. The substrates were 
then transferred to a different sputtering machine where a 
sputter etch of the target and substrate was performed 

O 

followed by the deposition of approximately 9,000 A of Si0 2 . 
That system's base pressure was typically 2 x 10“® torr. This 
layer was subsequently patterned with a photolithographic 
step in the base electrode geometry, as shown in Fig. 5.2. 
Utilizing the photoresist as a mask, the exposed Si0 2 was 
removed with a BHF etch and the photoresist subsequently 
removed with acetone. An ion milling process was then 
performed, using the dielectric pads as a mask to the 600 eV 
argon ion beam, removing the exposed Nb film. This step 
creates an edge in the Si0 2 / Nb structure on which the 
junction area will subsequently be defined. The substrate 
mount is cooled to 200-250 K during this step for two 
reasons. Without cooling the ion beam cross links the 
photoresist, making subsequent removal of the photoresist 
very difficult. It has also been shown that, without 
cooling, ion milling at these energies will adversely effect 
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side view top view 



Figure 5.2 Base electrode patterning with photoresist 


side view top view 



Figure 5.3 Gold contact pads on Nb 
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the superconducting quality of the Nb film. 

The next photolithographic step defined a window in the 
photoresist on the back half of the base electrode. The 
exposed Si02 was then removed with a BHF etch, forming a 
contact window to the Nb substrate. A 4,000 A thick film of 
gold was then RF sputter deposited onto the chip after a 
sputter etch cleaning of the exposed Nb surface. The regions 
of gold covering the photoresist were lifted off with 
acetone, leaving gold contact pads on the Nb surface (Fig. 
5.3) . 

The final photolithographic step involved definition of 
the photoresist counterelectrode liftoff stencil. This 
stencil primarily consisted of a thin open 'finger', 2-4 urn 
wide, overlapping the Si0 2 /Nb edge with an undercut in the 
profile of the photoresist as shown in Fig. 5-4. A 600 eV 
argon ion beam was then employed to clean the exposed Nb, 
using the photoresist stencil as an ion milling mask. The 
substrate mount was cooled to approximately 200 K during this 

o 

step. Following the removal of 100-200 A of surface oxide 
and Nb, such that the edge plane presents a clean Nb surface, 
the ion beam was switched off. 

An approximately 30% 0 2 , 70% Ar gas mixture was then 
obtained and the 600 volt ion beam re-established within a 
period less than one minute. This reactive ion beam 
oxidation step was then employed for approximately 30 seconds 
to form a thin Nb 2 05 (native oxide) tunnel barrier on the 
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exposed Nb. A PbBi counterelectrode was then immediately 
deposited by thermal evaporation, after a 90% rotation of the 
substrate holder away from the ion beam, toward the 
evaporation source. A pre-weighed 28% Bi and 72% Pb charge 
was evaporated to completion such that an approximately 4,000 
A thick film was obtained. This particular PbBi composition 
is generally referred to as e phase PbBi, and results in a 
single phase alloy film. The evaporation system was equipped 
with a shutter such that the PbBi can be outgassed and heated 
to a liquid state before counterelectrode deposition. A 
liftoff step was performed, after removal of the chip from 
the vacuum system, by soaking the chip in acetone, resulting 
in five or six completed devices. The processing was 
completed with the use of a protective covering layer of 
photoresist to prevent further oxidation and adsorbtion of 
moisture from the atmosphere. A more detailed discussion 
concerning device development and changes in fabrication 
technology from January 1982 to April 1984 follows. This 
material will be presented for the most part chronologically 
such that different aspects of the general fabrication 
technology will be examined more closely as its relevance 
becomes apparent. 

Initial results obtained with this technology were very 
promising. In Fig. 5.5, a typical I-V curve of a Nb/PbBi 
device (1*2,4 = 6.5) is shown. Though the general tunneling 
characteristics of an ideal SIS device are present, many 
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non-ideal effects are also exhibited. Of particular 
importance is the fairly rounded quasiparticle current jump. 
The I-V 'turn on' structure, expected to be an abrupt 'knee' 
atA^ + A 2 = 3.14 mV, is actually quite rounded for biases 
between 1 and 3 mV. However, the subgap current, for biases 
below 1 mV, is relatively low. 

These non-ideal effects indicated the existence of a 
distribution of gap energies. It was hypothesized that the 
superconducting quality on the edge of the Nb film had 
deteriorated. It is well known that value of the energy gap 
of Nb decreases for films less than 1,000 A thick. Utilizing 
an Si 02 ion milling mask for formation of the base electrode 
edge (as described above), the typical edge angle achieved 
was less than 30°, as shown in Pig. 5.6. Nb also has a 
rather short coherence length of -100 A. For 30° edge 
angle, there would exist a film edge length of approximately 
200 A, which would be supported by 100 A of Nb (Fig. 5.6). 
It was, therefore, thought that the tapered portion of the Nb 
edge was presenting a graded distribution of energy gap 
values, resulting in a rounded current jump below the 
expected sum potential. 

The choice of Si0 2 as a milling mask for the base 
electrode step had been made for a number of reasons. The 
technology for Si0 2 deposition and the patterning of such 
films by chemical etching had been previously developed in 
the Semiconductor Device Lab. Unfortunately, the ion milling 
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Figure 5.6 Tapered edge profile with edge angle of 30° 


Al a°3 

Nb 


Figure 5.7a AI 2 O 3 pad on top of the Nb film 



PR 


Nb 


Figure 5.7b Photoresist definition prior to base 

electrode formation with ion miling 
machine 



Figure 5.7c 


Resulting A^C^/Nb base electrode 
edge structure 
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rates of SiC>2 and Nb films are nearly identical, resulting in 
a small edge angle. In order to obtain a steeper edge angle, 
an alternative AI2O3 milling mask technology was developed. 
The milling rate of AI2O3 is more than three times less 
than the rate of Nb such that the use of AI2O3 in the edge 
milling step should result in a more desirable, steeper 
profile. 

This process was begun with the deposition of 2,500 A 
of A1 by RP sputtering in an argon plasma. A 
photolithographic step, which is nearly the inverse of the 
previous base electrode step, was then performed. An open 
window in the photoresist, of the same general shape as the 
previously solid photoresist base pads, was defined on top of 
the Al/Nb metal layer. The exposed A1 layer was then 
converted to A1 2 C>3 by an anodization process. The 
photoresist and remaining A1 were then removed, leaving an 

AI2O3 pad on top of the Nb film (Fig. 5.7a). Since AI2O3 is 
very difficult to pattern chemically, the length of this pad 
was designed to be shorter than the previous, full length 
base electrode structure. Another photolithographic step 
was, therefore, executed, defining solid photoresist pads 
which partially overlapped the Al 2 o 3 structures, as well as 
regions of the Nb film (Fig. 5.7b) to allow eventual 
electrical contact for device testing. An ion milling step, 
exactly the same as used in the Si0 2 technology, was employed 
to form the final A^C^/Nb base electrode edge structure 
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(Fig. 5.7c). 

In order to develop this new process, a number of 
technological problems had to be solved. Some of the more 
sensitive processing parameters will, therefore, be discussed 
below. 

The major difficulty encountered centered on the 
anodization of the exposed A1 film. The anodization was 
performed in a 0.05 M ammonium tartrate solution with the 
aluminum film as the anode and a strip of platinum for the 
cathode. The anodization cell was driven by a constant 
current supply which yielded a current density of 
approximately 500 yA/cm^. The anodization rate was 
calibrated at ~13 A/V. Unfortunately, the thickness of the 
resulting AI2O3 film is 1.7 times the thickness of the 
original fim. This differential resulted in poor adhesion 
of the photoresist to the growing AI2O3/AI step. As the 
thickness of the AI2O3 film was increased, the edges of the 
photoresist mask would pull away, resulting in a ragged AI2O3 
definition. 

Initial results for thermally deposited A1 films were 
very poor as less than 500 A of A1 2 0 3 could be formed before 
the photoresist would begin pulling away. If the anodization 
was continued, the entire layer of photoresist would 
eventually lift off. This problem was alleviated with three 
changes in the fabrication technology. First, the method of 
Al deposition was changed from thermal evaporation to RF 
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sputtering. If the substrate was thermally attached to the 
cooled sputtering table, the resulting small grained films 
were much harder and more durable than the evaporated films. 
Secondly, it was determined that if the traditional organic 
cleaning step (ethanol, tricloroethylene and methanol) 
performed prior to base electrode photolithography was 
omitted then the photoresist adhesion was greatly improved. 
The photolithographic step was, therefore, begun immediately 
following completion of the A1 sputtering. The third 
important step was the optimization of the photoresist 
postbake process prior to anodization. It was determined 
that liftoff was minimized if a 120 C postbake was used for 
20 minutes, the temperature of the postbake being the more 
crucial factor. 

Although the A1 2 0 3 technology successfully yielded a 
much steeper edge angle, the quality of the resulting I-V 
curves was not improved. In fact, the 1*2,4 rat i° for such 
devices (Fig. 5.8) had decreased from 6.5 to 5.3 with respect 
to the Si0 2 devices. Further work spent experimenting with 
the A1 2 0 3 technology resulted in devices with even lower R 2 ,4 
ratios. A return to the previous Si0 2 technology showed no 
improvement in device quality. The conclusion drawn was that 
the edge slope angle was not the dominant factor, i.e., there 
was some other fabrication step to which the device I-V 
characteristics were highly sensitive. This deterioration of 
device quality and the efforts involved toward pinpointing 
the source or sources of this degradation will be the subject 
of Chapter VI. 




CHAPTER VI 


A Search for the Mechanism of 
Device Quality Deterioration 

Despite attempts to identify and reverse the 
deterioration of device quality f the degradation continued 
until February 1983, at which point the I-V characteristics 
were completely ohmic, lacking all SIS structure. The 
efforts involved in tracking the causes of device 
deterioration and the subsequent alterations to equipment and 
fabrication procedure will be discussed in this chapter. 

Before the degradation search is considered, it would be 
helpful if Appendices A and B were examined. In Appendix A, 
the general fabrication procedure is explained in more detail 
than outlined in Chapter V. Specific temperatures, times and 
chemicals, etc. used in the processing steps are given for 
the reader to refer to when considering subsequent processing 
changes. A new Si0 2 base electrode defining technology was 
developed after the A1 2 0 3 method was optimized. This new 
Si0 2 method utilizes a photoresist/Si0 2 ion milling mask, as 
opposed to just the Si0 2 of the previous structure, to define 
the base electrode geometry. Since this technology was much 
simpler than the A1 2 0 3 method and resulted in fairly steep 
edge angles, the new Si0 2 method was utilized in most of the 
remaining SIS efforts. 

In Appendix B, an introduction to the theory of ion gun 
operation is given. Since the extremely sensitive device 
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fabrication steps of cleaning and oxidation involve the use 
of an ion gun technology, a knowledge of ion gun operation is 
imperative. A description of the ion gun and associated 
vacuum chamber used for much of this work is also included in 
Appendix B. 

In the search for the degradation mechanism one of the 
first areas pursued was the integrity of the superconducting 
materials. Traditionally, the characteristics of the 
superconducting transition of the material is used as a 
quality test for thin films [28]. A T c value below the 
transition temperature established for bulk material or a 
gradual superconducting transition is an indication of poor 
quality material. Both the Nb and PbBi showed sharp 
transitions at 9.1 K and 8.3 K respectively. The critical 
current densities for both materials were obtained from weak- 
link structures. The measured values for critical current 
density of 10 7 A/cm 2 correspond to bulk values, indicating 
quality materials. • An alternative quality parameter, often 
quoted for Nb, is the ratio of a film resistance at room 
temperature to the resistance measured at 10 K (commonly 
referred to as the resistance ratio). The measured value of 
~3 corresponded to values commonly reported for Nb [18] . 

Though these results were positive, these tests had not 
originated as a test of quality for SIS tunnel devices but 
rather only as an indicator of the superconducting quality of 


the material. 
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The possibility that a contaminant was being introduced 
in one of the processing steps was also fully examined. The 
extremely crucial step of ion cleaning and oxidation was 
first investigated. All chamber walls, in a direct line-of- 
sight with the ion beam, were covered with pure Nb foil. A 
Nb cover for the copper substrate holder was also installed 
in an attempt to reduce possible backsputter ing of impurities 
during milling. The ion milling gun, chamber and diffusion 
pump were disassembled and thoroughly cleaned. This cleaning 
process was composed of several steps: an acetone wash, a 

hot H 2 0 and alconox wash, several hot H 2 0 rinses, a TCE 
rinse, a methanol rinse and a 200°C oven bakeout. All 
O-rings were replaced and the diffusion pump oil changed. 

During the cleaning, copper deposits were found on the 
two alumina insulating rings adjacent to the filament of the 
ion gun. This copper was deposited by low energy sputtering 
of the two copper pieces supporting the filament cathode. 

This design flaw was corrected by mounting the filament via 
two electrical feedthrough insulators on a brass back plate 
which is in direct physical contact with the anode chamber. 
The cathode-anode potential difference is thereby established 
only at the filament and its holders, significantly reducing 
the sputtering of metallic impurities from the interior of 
the gun. The original and replacement design are shown in 
Fig. 6.1. 
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It was also discovered that the carbon grids, provided 
with the gun, were formed from an unpurified grade of 
graphite. These grids were replaced with pyrolitic graphite 
(obtained from B.F. Goodrich Company) containing less than 5 
ppm total impurities. Since the ion filament is constantly 
subjected to ion bombardment (~ 50 eV) , the filament was 
periodically replaced in an attempt to return the system to 
the initial operating conditions. Though it was felt that 
all of the above modifications to the ion gun were 
beneficial, no improvement in the device quality was seen. 

In an attempt to isolate the effects of the ion 
oxidation step, a number of devices were fabricated utilizing 
a thermal oxidation step to form the Nb 2 0 5 tunnel barrier. 

In this method the filament and the argon were turned off 
after the ion cleaning process and an oxygen pressure of 
4 x 10”4 torr was established in the chamber. This pressure was 
maintained for a period of 60 minutes while a thermal oxide was formed. 
The chamber was subsequently evacuated (pressure < 1 x 10"*® 
torr) and the substrate holder, which had been allowed to 
thermally drift from approximately 180 K to approximately 240 
K, was recooled to a temperature T ~ 180 K. The standard 
process of evaporation was then resumed. A typical I-V curve 
resulting from this method is shown in Fig. 6-2. 

Unfortunately, this test did not alter the continuing trend 
of device quality degradation. 



Figure 6.2 



Typical I-V curve obtained with thermal 
oxidation technology 



Figure 6.3 


Typical I-V curve obtained with extended 
cleaning step (40 min.) 



87 


The peculiar nature of this problem was a very slow rate 
of deterioration. The quality of devices fabricated on chips 
on successive months deteriorated slightly but steadily. 

The impetus for a number of changes to the fabrication 
schedule was centered on the hypothesis that some sensitive 
portion of the various fabrication processes was slowly 
becoming less effective. For example, the duration and 
temperature of the postbake used in the photolithographic 
base electrode step was re-examined. Devices were fabricated 
with postbake parameters changed by varying degrees from the 
standard method. 

Another consideration was the ion milling cleaning step, 
prior to oxidation. Cleaning times of various durations were 
examined. It was found that cleaning times less than 
approximately 5 minutes resulted in a linear I-V structure, 
however, cleaning durations 2 to 3 times the nominal period 
of 15 minutes yielded no improvements (Fig. 6-3). 

Lower ion beam supply voltages were also used in an 
attempt to minimize possible surface damage from the argon 
ions. The edge structure was still formed with a 600 eV 
beam, however, beam supply votlages of 300, 400 and 500 
volts, for the cleaning and oxidation step, were utilized. 

In Figure 6-4, the I-V curve resulting from such a run with 
V BS “ volts is shown. There is no trace of SIS 
structure, the I-V characteristics were simply ohmic. At 
this juncture the task of tracing the source of degradation 
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Figure 6.4 


Typical I-V curve obtained with low voltage 
cleanand oxidation (350 volts) 
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became even more difficult. With no SIS I-V structure/ only 
"positive" effects on device quality could be measured, 
whereas previously both "positive" and "negative" effects 
could be seen. Unfortunately, during much of the work after 
this point, there was little evidence that various 
alterations in the fabrication technology were having any 
effect, positive or negative. 

Numerous efforts were also focused on the possibility 
that an undesirable contaminant had been introduced through 
one of the reagents or materials used in the fabrication 
steps. All processing solvents, etches and photolithograhic 
chemicals were replaced. Different ovens were used in the 
processing and all the beakers were replaced. A method in 
which the Pb and Bi pellets were alloyed in a pure carbon boat 
was developed. The pellets were alloyed in a forming gas 
atmosphere (95% N2/5% H 2 ) for 5 minutes at which point the 
charge was well mixed and shiny. The charge was then removed 
and immediately placed in a beaker of methanol to prevent 
formation of a surface oxide. The charge was 
subsequently set in the evaporation boat and the chamber was 
quickly roughed out and opened to the diffusion pump. 

As a further check on the quality of the superconductive 
materials, an x-ray analysis was performed on an RF device. 

No impurities were detected except Ar in the sputtered Si0 2 , 
a common occurrence for films sputtered in Argon plasmas. 
Unfortunately, this technique is not sensitive to elements 



with atomic weights below 24 - most notably carbon and 
oxygen. 
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One rather interesting hypothesis centered on the 
possibility that a crucially needed "contaminant" present 
initially in some processing step, had been slowly removed or 
depleted. In particular, a 'contaminant' present in the ion 
milling gun or chamber might have been removed through 
repeated use. The system might have become too clean! 

An analogous case involving the fabrication of large 
area SIS devices (where the presence of an unintended 
impurity was necessary for optimal results) had been reported 
by Raider et al. [29,30] Their work involved fabrication of 
Nb/Pb alloy overlap devices, utilizing an argon RF plasma 
etching of the Nb film prior to oxidation. An oil diffusion 
pump was used to evacuate their system to ~10“® torr while 
the cleaning and oxidation steps were performed at a pressure 
of ~3 utorr. The conditions required to obtain optimal 
tunnel barriers with this system were quite unexpected. 

Tunnel barriers with the most nearly abrupt Nb/ltt^C^ 
interface yielded I-V characteristics with the highest subgap 
current. However, barriers with a 25-30 thick NbCO 
intermediate region between the Nb and the native (Nb 2 C>5) 
oxide [31] resulted in optimal devices. The NbCO formation 
was traced by Raider to the plasma cleaning step where the 
thickness of this transition region was found to be directly 
related to the applied RF voltage. Utilizing RF voltages of 
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-1,100 V, oxygen and carbon impurities, which were 
unintentionally present in the chamber, were implanted into 
the Nb. Subsequent plasma oxidation of this surface 
transformed a portion of the NbCO region into the ^205 
layer. Devices without the NbCO layer were found to have 
regions of NbO between the Nb and the Nb205 tunnel barrier. 

The interpretation given for these results [32] was that 
the NbCO layer shielded the Nb from direct oxidation, 
preventing formation of NbO - a conductive suboxide of Nb. 

The NbCO was thought to act as a transition region between 
the oxygen sensitive Nb and the desired Nb2C>5 barrier. 

Returning to the degradation problem, the possibility existed 
that a carbon impurity, present initially in the ion gun 
chamber and required for formation of a NbCO transition 
layer, had been slowly cleaned away. A lecture bottle of 
methane gas, to be used as a source of carbon impurity, was 
incorporated in the gas manifold. Several different 
schedules for introducing the CH4 during the cleaning and 
oxidation steps were developed. In one such schedule the 
cleaning step was performed with a gas mixture of 10% CH4 and 
90% Argon. In another schedule the CH 4 was admitted for only 
the last minute of the fifteen minute cleaning step. These 
and other methods, however, failed to change the resulting 
ohmic I-V characteristics of the SIS devices. 

The growing concern about the background gas composition 
led to a simple test of vacuum integrity. A chip was 
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processed following the usual procedure until the ion 
cleaning step was completed. At this point one of two 
different procedures would be followed. The first method 
involved the omission of the oxidation step by immediately 
turning the substrate into a flux of PbBi. The second method 
involved the omission of the oxidation step, followed by a 
one minute pause before evaporation of the PbBi (which formed 
the counterelectrode). Critical currents were then measured 
for the resulting devices. Critical current density values 
resulting from the first method corresponded to bulk values. 
Results from the second method, however, were two order of 
magnitudes lower than bulk values. During the one minute 

o o 

pause approximately 8 A - 12 A of background gas had formed 
on the surface. 

It was not known whether this condition had always 
existed in the system or if it was symptomatic of the 
degradation problem. A leak detector was not conveniently 
available to check the integrity of the vacuum system, 
however, an unconventional method was utilized to initially 
check the system. The vacuum system was evacuated to 9 x 10“ 
torr. Helium gas was then directed toward the various 
seals of the vacuum chamber. When the gas was aimed at the 
ion gun, the pressure indicated by the ionization gauge 
decreased to 6 x 10" 7 torr. It was apparent that the He gas 
was entering the chamber through a leak in the ion gun and 
effecting the ionization gauge operation. It was determined 
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that the O-ring seals, associated with each insulator ring, 
leaked (see Appendix B). 

The sites of the leaks were in a very sensitive region 
since entering gas impurities would diffuse directly into the 
ion beam. It was discovered that the o-ring grooves on the 
magnet, two grid rings, and the two neutralizer rings had 
been improperty machined. After a replacement of the o-rings 
and a machining of the grooves failed to improve the vacuum 
seal, it was decided to construct a new vacuum chamber in 
which to entirely enclose the existing ion gun. A chamber 
was designed (Fig. 6.5 and 6.6) which had a larger clearance 
between the beam path and chamber walls and which utilized 
only nine o-rings in comparison to twenty-five for the 
previous chamber. The ion gun was mounted directly onto a 
water cooled brass back plate which was isolated from the 
chamber ground by a plexiglass ring. The o-rings associated 
with the insulating rings of the ion gun, providing no 
further purpose, were removed. 

Nb/PbBi junctions were subsequently fabricated with the 
new chamber. The resulting I-V device curve, though of 
rather poor quality with an excessive subgap current, 
possessed a weak SIS tunneling structure with a turn-on 
voltage at the sum of the gap energies (Fig. 6.7). The 
enclosure of the ion gun within a separate vacuum chamber 
reversed the pattern of device quality degradation for the 
first time in two years. However, it appeared that enclosing 







Figure 6.6 


New vacuum chamber design 
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the ion gun only slightly diminished the effect of some 
contamination source which was effecting the quality of the 
tunneling barrier. The integrity of the insulator rings as a 
vacuum material was immediately suspected. Five plate glass 
rings were fabricated by slowly grinding the pattern in 
pieces of plate glass with two concentric brass tubes. 

With the new glass rings installed, the operational 
pressure range increased from a 1 x 10”^ - 4 x 10 - ^ torr 
interval to a 5 x 10“^ - 4 x 10”^ torr pressure 
range. Nb/PbBi junctions were subsequently fabricated? the 
resulting I-V curves of two junctions, representative of the 
chip, are shown in Fig. 6-8. A dramatic improvement in 
device quality was evident? the subgap current was much 
smaller and the turn-on voltage knee was sharper. The 1*2,4 
ratios of these two devices were 3.8 and 4.7. Though these 
ratios were less than the highest ratios obtained (6.5) 
before the degradation problems occurred, these results were 
very encouraging. The old alumina ion gun insulating rings 
were thought to have been poorly fixed, such that a porous 
interior material was exposed after a year of ion milling. 

These rings, which collar the plasma and the existing ion 
beam, were outgassing (acting as a virtual leak) directly 
into the plasma and ion beam, contaminating the tunnel barrier. 

Two further chips were subsequently fabricated in order 
to verify the results of chip #37. The I-V curve of the 
highest quality junction is shown in Fig. 6-9. The 1 * 2,4 





100 


ratio was equal to 9, at that time this was the largest ratio 
ever obtained for a SIS device at the University of Virginia. 
Unfortunately, devices of only fair quality, characterized by 
r 2,4 rat ^ os as l° w as 2.4 , were fabricated on the same chips. 
It was suspected that the photoresist counterelectrode might 
have been backsputtered into the junction area during the 
cleaning and oxidation steps, causing nonuniform results. A 
device structure, utilizing a Si0 2 window (collared by a 
photoresist counterelectrode) to define the device area, was 
subsequently developed (Fig. 6.10). Though two chips were 
fabricated, no improvement of electrical uniformity across 
the chip was found. 

The batch of Nb, which was being utilized at that time, 
was examined and the value of the resistance ratio was found 

^ R roon/ R 10K) t0 have decreased from a previous value of three 
to a significantly lower value of two. These films were 
thought to be of too low a quality for SIS tunnel junction 
fabrication such that an alternative source of Nb films was 
sought. 

At a meeting with Dr. d'Addario of NRAO, it was 
suggested that James Beall [33] at NBS, who was working on Nb 
edge junctions in Boulder, Colorado, might be a potential 
source of Nb films. Upon contact, Mr. Beall very kindly sent 
four wafers of Si on which 2,500 A of Nb had been DC 
magnetron sputtered. The resistance ratio value for these 
films was 3.5. The electrical quality of junctions 
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fabricated on a chip of this material was highly uniform with 
1*2,4 rat i° s of approximately 5.6. A typical I-V curve is 
shown in Fig. 6.11 and numerical values for R 30 O' R 2' R 4 an< ^ 
R 2 f 4 are tabulated in Table 6-1. A typical I-V curve 
obtained by Mr. Beall on this same Nb is shown in Fig. 6.12. 

A comparison of Fig. 6.11 and 6.12 indicates that the 
fabrication technology used in work described herein is 
superior. Unfortunately, the Si substrate material has very 
poor rf characteristics which prevented efficient coupling 
into such devices when mounted in mixer substrates. 

Subsequent Nb films, deposited on quartz substrates at NBS, 
proved to be of inferior quality, yielding poor, nonuniform 
I-V characteristics. 

At this point, a great deal of confidence existed in our 
fabrication technology and device design; however, a source 
of high quality Nb films did not appear to be forthcoming. A 
switch from Nb to NbCN as a base electrode material was 
subsequently made. NbCN has a superconducting transition 
temperature of ~ 16 K, a large energy gap of '-2.4 eV and 
presents a very physically robust film. Of extreme 
importance is the reported characteristic of a high quality 
native oxide, with the absence of conductive suboxides. [34] 
Dr. Cukauskas of NRL [35] was contacted and he very nicely 
agreed to supply test NbCN films RF sputtered onto quartz 
substrates. This work is presently continuing toward 
optimization of NbCN/Pb alloy devices. 














Figure 6.11 


Typical I-V curve obtained with NBS Nb 



Figure 6.12 


Typical I-V curve obtained by Mr. Beall 
on NBS Nb 



CHAPTER VII 


Investigation of NbCN 

The switch to a NbCN base electrode did not require a 
significant alteration of the fabrication technology 
developed for Nb based junctions. The milling rates of the 
two materials are approximately equivalent, neither material 
is etched by BHF acid and both materials retain their 
superconducting properties when deposited as films of 
thickness greater than 1,000 A. The anodically formed oxide 
of NbCN, however, is physically fragile and, therefore, 
cannot be utilized in a device structure. 

Initial development with the ion gun technology yielded 
junctions which were electrically uniform across the chip, 
possessing low subgap current for voltages approximately 
two-thirds below the gap voltage (Fig. 7.1). Though the 
uniformity of results and low subgap current was very 
encouraging, the quasiparticle current jump was very 
rounded, indicating a damaged surface layer of NbCN. It was 
thought that the formation of the edge structure by ion 
milling was significantly damaging the superconductor 
lattice, resulting in a distribution of gap energies across 
the junction area. The ion cleaning and oxidation of this 
junction was performed at 450 volts. Junctions were then 
fabricated utilizing a 350 volt ion beam, resulting in 
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devices with improved, sharper tunneling characteristics at 
the gap voltages. 

Since the ion gun used in this work was designed to 
operate at 600 eV, a number of modifications were made to 
increase the low beam currents obtained for voltages below 
300 eV. Junctions were then fabricated, utilizing a 200 eV 
ion beam, with measured electrical characteristics which 
were superior to all previously fabricated devices at this 
lab (Fig. 7.2). Unfortunately, the efficiency of the ion 
gun used in this work deteriorated rapidly below 200 eV. 
Junctions have been fabricated with an ion beam voltage of 
125 volts, however the resulting electrical characteristics 
were poor and non-uniform across the chip. Occasionally, a 
single device with high quality electrical characteristics 
was fabricated (Fig. 7.3) on a chip with generally poor 
electrical characteristics. It is therefore believed that 
further improvement of device electrical characteristics 
could be obtained if a high current density, low voltage 
ion beam was utilized. 

After numerous attempts toward improving the 
performance of the ion gun had failed, a new 3 cm ion 
source (Model I) was obtained from Commonwealth Scientific 
Corporation. The efficiency of this ion source, when 
operated with a beam voltage greater than 100 eV, was 
excellent. Devices fabricated with this ion source, 
however, were of poor quality, even when the fabrication was 




Figure 7.2 



Figure 7.3 


performed with a 400 eV ion beam voltage. The base pressure 
of the chamber was 5 x 10~ 8 torr, however after 20 minutes 
of ion source operation the pressure would rise to 1 x 10 -8 
torr. Again the insulators were immediately suspected, 
however due to the geometry of the ion gun it was impossible 
to fabricate glass replacement pieces. Tungsten, thoriated 
tungsten, tantalum and Nb filaments were all examined in an 
attempt to find a more efficient emission source. Different 
ion source operation schedules have also been examined. For 
one such run, the ion gun was operated at a much higher 
discharge and filament current than normal, with the sample 
turned away from the beam. The ion source was then turned 
off and the pressure allowed to return slowly to a 10 -8 torr 
range. This process was repeated five times in an effort to 
rid the insulators of contaminants. The junctions which 
were subsequently fabricated possessed typically poor 
electrical characteristics. 

In an effort to provide an efficient low voltage ion 
source, the vacuum chamber was modified to permit plasma 
etching (Fig. 7.4). The plasma technology has been utilized 
successfully for large area planar NbCN junctions (Hikita). 
It should also permit a low voltage (hence low damage) 
cleaning and oxidation of the NbCN edge surface. Initial 
tests of the plasma system indicate that a plasma can easily 
be obtained for discharge potentials above 25 volts and gas 
pressures above 1 pm. Devices will be fabricated after the 
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etch rate of the system as a function of discharge voltage 
is determined. It is expected that the most difficult 
parameters to determine will be the oxygen partial pressure 
and the oxidation time. 


CHAPTER VIII 


Conclusion and Recommendation for Future Research 

During the course of this work, a number of procedures 
were developed for fabricating Nb/Pb alloy and NbCN/Pb allow 
devices. The quality of the tunnel barriers resulting from 
the oxidation of Nb has been found to be extremely sensitive 
to contamination. A great deal of effort was expended 
toward solving a critical degradation of I-V 
characteristics, due to barrier contamination. It was 
discovered that the electrical quality of the device was 
even sensitive to the virtual leaking of gases from porous 
ion gun insulating rings. It was also shown that the device 
I-V characteristics are extremely sensitive to the quality 
of the Nb material. 

An investigation of the use of NbCN base electrode 
films is presently underway. The reported [36] absence of 
conductive suboxide in the native oxide of NbCN is very 
encouraging. However, the superconductive quality of the 
NbCN surface has been found in this work (and others [37]) 
to be very sensitive to ion bombardment and subsequent 
surface damage. The use of a plasma etching technolology is 
presently being explored. This will hopefully permit a low 
voltage cleaning and oxidation of the NbCN edge surface. 

The introduction of CF 4 during the cleaning process in order 



to clean "chemically" instead of "physically" (perhaps 
resulting in less surface damage) is also being considered. 

The use of an artificial barrier should also be 
investigated. Excellent results have been obtained by 
Gurvitch et. al. [38] using an AI 2 O 3 barrier with 
60 m 2 - 400 ym 2 junction area Nb/A^C^/PbBi and all Nb 
structures. The subgap current of these devices have also 
been fabricated [39] which have an unprecedented I-V 
quality, comparable to the I-V quality of the best Pb based 
junctions. However, successful micron sized junctions have 
yet to be fabricated. 

The plasma, CF^ and AI 2 O 3 technologies are presently 
being developed at this lab. The initial results of the 
plasma technology are promising, with ** 2,4 rat i° s greater 
than 20 measured for 2.5 ym 2 devices. Initial results from 
the AI 2 O 3 technology are expected before January. It is 
felt that all three technologies hold great potential and 
are worth pursuing. We are presently investigating 
alternative avenues for future funding. 



Appendix A 


The general procedure used for fabricating Nb/PbBi 
tunnel junctions of ~1 ym^ area is summarized in this 
appendix. While a particular sequence of specific steps is 
given, variations on each were performed in the course of 
development of this fabrication technology. 
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I. Mi BE Deposition 

A. Clean quartz substrates 

1. Boil and degrease in ethanol, 
trichloroethylene, and methanol ("organic 
clean") 

2. Dry bake 120 °C for 20 min. 

B. RF sputter deposition of 2,000 A of Nb 

1. System base pressure of 2 x 10~ 7 torr 

2. Deposition rate of 40 A/sec., 10 urn Ar 

3 . H 2 O cooled substrate table 

II. £Li &2 EE De position 

A. Clean chip 

1. organic clean 

2. Dry bake 80 °C for 20 min. 

B. RF presputter of 3,000 A Si0 2 in "empty" chamber 

C. RF sputter deposition of 3,000 A of Si0 2 

1. System pressure of 1 x o 10 -6 

2. Deposition rate of 50 A/sec., 15 ym Ar 

3. H 2 0 cooled substrate table 

4. 500 volts (target), 50 volts (substrate) 
hi. Bas.e, Electrode PhPtoli.tho g x.a p hy 

A. If this photolithographic step follows immediately 
after Si0 2 deposition, then the following cleaning 
step is omitted. 

B. Clean chip 

1. Organic clean 

2. Dry bake 80°C for 20 min. 

C. Spin AZ1350J photoresist 

1. 3,500 RPM 30 sec. (approximately 2 ym) 

D. Treatment prior to illumination 

1. Prebake 80°C for 25 min. 

2. Allow to cool at room temperature (5 min) 
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E. Exposure with base electrode mask 

1. Relative UV intensity 49 (Preco Aliguer) 

2. 25 sec. 

F. Develop photoresist 

1. AZ developer/de-ionized water (DI) 1:1 

2. Develop by "inspection" (The photoresist is 
developed until it appears to have stopped 
being removed ( ~ 20 sec.). Development is then 
continued a few seconds further.) 

3. Three rinses in dionized water (DI) 

4. Postbake 120°C for 20 min. 

iv. inn Mill nf Rase Electrodes 

A. Chip mounted with Apiezon-L grease to Nb covered 
substrate holder 

B. Operate ion gun for 30 min. with substrate turned 
away from beam 


1 . 

V BS = 600 volts 


2. 

V A = -50 volts 


3. 

V dis = volts 


4. 

I fil = 14 am P s 


5. 

I b = 50-70 yA/cm^ 


6. 

I dis = 0.1 - 0.3 amps 


7. 

System pressure <. 1 x 10“ 6 torr 


8. 

Cool substrate holder to 180 - 230 K 

of ln 2 

by means 

Ion 

mill base electrode 


1 . 

Continue milling 15 minutes past the 

time at which 


the exposed Si0 2 and Nb are removed 
substrate becomes visible) 

(quartz 

2. 

Warm chip to room temperature by forcing N 2 9 as 
through LN 2 feedthroughs 
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D. Treatment after ion milling 

1. Rinse twice in boiling xylene to remove 
Apiezon-L grease 

2. Remove photoresist with aceton 

3. Rinse with methanol 

V. Contact Window Format i-Sn 

A. Clean chip 

1. Organic clean 

2. Dry bake 80°C for 20 min. 

B. Spin AZ1350J photoresist 

1. 4,000 RPM 30 sec. 

C. Treatment prior to illumination 

1. Prebake 80 C for 20 min. 

2. Allow to cool at room temperature (5 min.) 

D. Exposure with contact window mask 

1. Relative UV intensity 49 (Preco) 

2. 25 sec. 

E. Develop photoresist 

1. AZ/DI 1:1 

2. Develop by "inspection" 

3. Three rinses in DI water 

4. Postbake 120°C for 20 min. 

F. Open contact windows in Si02 with BHF 

1. BHF 3-1/2 - 4-1/2 minutes (60 ml HjO, 40 ml 
NH 4 F, 9 ml HF) 

2. Two rinses in DI water 

VI. Gold Contact Deposl-tian 

A. Affix chip to H 2 O cooled substrate platform with 
Apiezon-L grease. 

B. RF sputter of Au with shutter covering chip 

1. System Q pressure of 1 x 10 -5 torr 

2. 5,000 A Au 
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C. Continue pumping until system pressure of 
~2 x 10"® is obtained 

D. RF sputter of Au with shutter covering chip 

1. 15 ym $r 

2. 3,000 A Au 

3. 700 volts (target), 75 volts (substrate) 

E. RF sputter etch of chip 

1. 15 nm Ar 

2. -100 A removed 

3. 500 volts (substrate), 50 volts (target) 

F. RF sputter of Au with shutter covering chip 

1. 15 ym Ar 

2. 1,000 A Au 

3. 700 volts (target), 75 volts (substrate) 

G. RF sputter deposition of Au on chip 

1. 15 ym Ar 

2. 6,000 A of Au 

3. 700 volts (target), 75 volts (substrate) 

H. Treatment after Au deposition 

1. Rinse twice in boiling xylene to remove 
Apiezon-L grease 

2. Liftoff gold by soaking in acetone 

3. Remove with methanol 

VII. Counter Electrode Photolithogra phy. 

A. Clean chip 

1. Organic clean 

2. Dry bake 80* C for 20 min. 

B. Spin AZ1350J photoresist 

1. 6,000 RPM, 30 sec (approximately 1.5 ym) 

C. Chlorobenzene soak 

1. 15 min. soak 

2. Blow dry with N2 gas and let sit for 5 mins. 
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D. Exposure with counter electrode mask 

1. Relative UV intensity 49 (Preco) 

2. 30 sec. 

E. Develop photoresist 

1. AZ/DI 1:1 

2. Develop for a period equal to twice the 
observed optical development time ( 90 sec.) 

3. Three rinses in DI water 

VIII. Xon Mill Clean./ Oxidation and PbBi Evaporation 

A. Load Ta evaporation boat 

1. 0.540 grams of Pb, 0.210 grams of Bi 

B. Mount chip with Apiezon-L grease to Nb covered 
substrate holder 

C. Operate ion gun for 30 min. with substrate turned 
away from beam 

1. V BS = 600 volts 

2. V A = -50 volts 

3. v di s = 50 volts 

4. = 14 amps 

5. 1^ = 50-75 yA/cm^ 

6. I< 3 i s = 0.1-0. 3 amps 

7. System pressure < 8 x 10“ 7 torr 

8. Argon pressure 1 x 10“ 4 torr 

9. Cool substrate holder to 180-210 K 

D. Ion clean 

1. 45° incidence 

2. clean for 15 min. 

3. Increase evaporation current to 40 amps 

4. Continue cleaning for 20 secs. 
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E. Oxidation and evaporation (this follows immediately 

after the cleaning step) 

1. Admit O 2 to total pressure of 1.3 x lO""^ torr 

2. I b = 20-25 yA/cm 2 3 4 

3. I fii = 15 amps 

4. Twenty-five sec. later, increase evaporation 
current to 60 amps 

5. Ten sec. later, turn gun off and increase 
evaporation current to 90 amps 

6 . Open evaporation shutter and turn substrate 
into evaporant 

7. Eyaporate to completion 

8. Warm chip to room temperature by forcing N 2 gas 
through LN 2 feedthroughs 

F. Clean chip and PbBi liftoff 

1. Rinse twice in room temperature xylene to 
remove Apiezon-L grease 

2. Liftoff gold by soaking in acetone (light 
ultrasonic agitation necessary) 

3. Blow dry with N 2 gas 

4. Spin protective layer of AZ1350J photoresist 


Appendix B 


The choice of an ion gun technology instead of a more 
conventional RF plasma technology was made for several 
reasons. With RF sputtering, the discharge is initiated with 
the application of an electric field. Excited electrons 
collide with the chosen working gas atoms resulting in an 
ionization of these atoms; once the plasma is established, it 
is maintained by secondary electron emission. Also, with RF 
sputtering, the substrates and the target are located 
essentially in the plasma. 

Ion sources, however, utilize thermionic emission to 
sustain the plasma. This permits operation of ion sources at 
pressures (< 10~ 4 torr) two orders of magnitude lower than 
conventional sputtering sources. Also, with ion sources, the 
plasma region is physically separate from the material to be 
milled. Because of these two characteristics, the ion source 
technology should be a cleaner process. 

In this Appendix, the fundamentals of ion milling 
operation will be outlined and the specific system utilized 
in this work will be considered. 

Section 1 

A typical ion source with power supplies is shown in 
Fig. B.l. The working gas of the system (typically Ar) 
enters the discharge chamber through a port in the back of the 


Filament 

supply 


+ Magnet Supply 



Figure B.l 


Typical ion source schematic 
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gun. Electrons from the filament are thermionically emitted 
and are accelerated (by a relativly low "discharge" 
potential, V^ s ~ 20-70 volts) to the anode chamber walls. 

Some of these energetic electrons from the cathode collide 
with the atoms of the working gas producing a plasma of 
electrons and ions. This ionization chamber is maintained at 
a potential V bs above ground. Some of these ions pass 
through the screen grid and are accelerated through the 
aligned holes of the negatively biased accelerator grid. The 
ion beam results from the combination of each of these 
individual beamlets. The plasma potential is approximately 
equal to the anode potential which is set by the beam power 
supply (V bs ). With a target connected electrically to 
ground, the energy of the ions impinging upon the chip is 
given by the beam power supply potential V bs » 

Since ion sources typically operate at such low 
pressures (~10“ 4 torr), the probability of an electron-source 
atom collision in the discharge chamber is typically 
increased by the application of an axial magnetic field 
between the anode and cathode. The much heavier ions are 
relatively uneffected. The discharge voltage (V^g) is 
typically chosen so as to produce only singly ionized atoms. 
Since a doubly charged ion will obtain twice the energy of a 
singly charged ion, it will also produce approximately twice 
the surface damage. Singly charged ions also strike 
the cathode with an energy equal to ev ds . 


Doubly charged 
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ions can, therefore, significantly increase the resulting 
sputtering contamination since the sputtering rates of most 
materials increase dramatically with ion energy above typical 
discharge potentials. 

The ion beam is extracted from the discharge chamber by 
a set of grids (the "ion optics"). Typically, the ion optics 
consists of two closely spaced graphite grids, though single 
and multigrid ion optic systems can be utilized for special 
applications. [40] The screen grid is not biased externally 
but allowed to float electrically. Its potential has been 
found to approximately equal the anode potential. The 
accelerator (or extractor) grid is commonly biased at a 
negative potential, V a , while the target is electrically 
connected to ground. The variation of potential in a typical 
dual grid system is shown in Fig. B.2. [41] 

The ion current is controlled by V T , the sum of the 
magnitudes of V a and V 5s (V T = V a + V bs ). The accelerating 
grid is negatively biased in order to obtain a larger ion 
current and to prevent backstreaming of electrons into the 
ion gun. The later reason is of particular importance if a 
neutralizer filament is positioned after the dual grids. In 
such a case, the filament is utilized to supply electrons to 
create a charge neutral beam such that the "neutralized" beam 
has equal densities of electrons and ions. The accelerating 
grid presents a negative barrier to these neutralizing and 
other electrons. Increasing | V Q | , however, will also 
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increase the beam divergence as well as the ion current 
striking the accelerator grid. The sputtering yield of this 
grid and resulting contamination of the target will also 
increase. The accelerating voltage should, therefore, be 
adjusted to avoid both backstreaming of electrons into the 
ion gun and excessive sputtering of the extractor grid. 

The ion current is also strongly controlled by the grid 
separation and geometry. Smaller grid spacing results in a 
larger electric field between the two grids (for a given V T ) 
and hence a larger ion current. Typical grid spacings range 
from 20-200 mils. [42] 

Section XI 

A schematic of the vacuum system, associated with 
the ion source utilized in this work, is shown in Fig. B. 3. The 
system is pumped by a diffusion pump through an optically 
dense liquid nitrogen cold trap. An ion gun, obtained from 
the Commonwealth Scientific Corporation, is mounted to the 
side of the aluminum chamber. This chamber is equipped with 
a LN2 cooled rotatable substrate holder, evaporation boat and 
evaporation shutter. The beam current density is measured by 
a 1 cm^ Nb covered current probe positioned approximately 
1/4" to the side and 1/2" behind the substrate holder. A 
quartz window is positioned on the top of the chamber to 
permit viewing of the substrate while milling. The working 
gases for this system are research grade argon and 
oxygen (Ar, 99.995%, 0 2 r 99.999%). All valves and tubes are 



Figure B.3 


Schematic of the vacuum system associated 
with the ion source in this work 
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stainless steel except a final 12" length of teflon tubing 
which is required for electrical isolation of the anode 
chamber. Fine metering valves are utilized to regulate gas 
flow into the anode chamber. 

A schematic of the ion gun is shown in Fig. Bl. An 
aluminum cyclinder serves as both the anode and a holder for 
the wire turns of the electromagnetic. The anode, screen 
grid, accelerating grid and retaining piece were each 
electrically separated by alumina insulating rings while 
O-ring seals exist between all of the components. Carbon 
grids with a 60 mil separation and 3 cm beam width are 
utilized. Both tungsten and thoriated tungsten filaments 
have been used with similar results. Typical electrical 
operating conditions for this ion gun are given in Table B-l. 
Subsequent changes made to this ion gun configuration are 
described in Chapter V. 
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System pressure 

1 X 

10“® torr 

Argon pressure 

1 X 

10“^ torr 

V BS 

600 

volts 

V A 

-50 

volts 

V dis 

50 

volts 

1 f i 1 

14 

amps 

I dis 

.2 

amps 

x b 

60 

IjA/cm^ 

■^sub 

200 

K 


Table B-l 

Typical operating conditions of the ion source system 
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